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The possible role of a first order QCD phase transition at nonvanishing quark chemical potential 
and temperature for cold neutron stars and for supernovae is delineated. For cold neutron stars, 
we use the NJL model with nonvanishing color superconducting pairing gaps, which describes 
the phase transition to the 2SC and the CFL quark matter phases at high baryon densities. We 
demonstrate that these two phase transitions can both be present in the core of neutron stars and 
that they lead to the appearance of a third family of solution for compact stars. In particular, a core 
of CFL quark matter can be present in stable compact star configurations when slightly adjusting 
the vacuum pressure to the onset of the chiral phase transition from the hadronic model to the NJL 
model. We show that a strong first order phase transition can have strong impact on the dynamics 
of core collapse supernovae. If the QCD phase transition sets in shortly after the first bounce, a 
second outgoing shock wave can be generated which leads to an explosion. The presence of the 
QCD phase transition can be read off from the neutrino and antineutrino signal of the supernova. 
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It has been realized in the last years, that the QCD phase diagram might exhibit a rich structure 
at high baryon densities, be it in the form of color superconducting phases ^ or the quarkyonic 
phase While recent lattice gauge simulations indicate that the QCD phase transition at van- 
ishing quarkchemical potential is most likely a crossover, there might exist a first order phase 
transition lines at nonvanishing quarkchemical potentials. This region of the QCD phase diagram 
entails properties of strongly interacting matter which are found in the core of neutron stars and 
core-collapse supernovae and will be probed by heavy-ion collisions with the CBM detector at the 
Facility for Antiproton and Ion Research (FAIR) at GSI Darmstadt. 

Neutron stars are bom in core-collapse supernovae as so called proto-neutron stars. The tem- 
peratures reached in those supernovae and in proto-neutron stars are up to 50 MeV with baryon 
densities well above normal nuclear matter density. Similar conditions are encountered in simu- 
lations of neutron star mergers. The masses of rotation-powered neutron stars, pulsars, have been 
quite accurately determined. More than 1700 pulsars are presently known, the best determined 
mass is the one of the Hulse-Taylor pulsar, M = (1.4414 ± O.OOO2)M0 [^], the smallest known 
mass is M = (1.18 ± O.O2)M0 for the pulsar J1756-2251 [§]. The most rehable lower limit for 
neutron star masses published in the literature is the one of the Hulse-Taylor pulsar (with only one 
noticeable exception [^, ^]). Note, that the mass of the pulsar J0751-I-1807 was corrected from 
M = 2.1 ±0.2Me down to M = 1.14 - 1.40Me The high masses and radii inferred from the 
X-ray burster EXO 0748-676 in an analysis done in [^] are not model independent, a multiwave- 
length analysis concludes that a mass of l.SSM© is more compatible with the data |^. In any case, 
high masses and radii do not exclude the possibility of having quark matter in the core of neutron 
stars [0]. 

In the following we will loosely denote the high-density matter as quark matter, although the 
QCD phase transition at high baryon densities is due to chiral symmetry breaking and not due to 
deconfinement. Let us first discuss the possible role of the QCD phase transition and the stability 
of compact stars in a toy model for quark matter with an equation of state of the form p = a - e 
with a constant a = 1/3 and a given energy density jump (see [11]). For the hadronic side, we 
use a relativistic mean-field model fitted to the properties of nuclear matter (here set GM3). If the 
phase transition occurs close to the maximum mass, always unstable solutions appear for the hybrid 
star with a quark matter core. The mass-radius curve changes its slope as soon as quark matter is 
present. For an onset of the phase transition at moderate densities, the presence of quark matter 
leads to stable configurations, the slope of the mass-radius curve does not change its sign [|T^, 11]. 

Color-superconducting quark matter can be described by the NIL model which includes both, 
dynamical quark masses via quark condensates and the color-superconducting gaps A for the three 
flavor case (see [ ll3| , 14] for astrophysically relevant calculations). The parameters of the model 
are the cutoff, the scalar and the vector coupling constants Gs, Gy, the diquark coupling Go, and 
the 't Hooft term coupling K. They are fixed to known hadron masses, the pion decay constant, 
which leave two free parameters. Go and Gy. In addition, the total pressure is usually fixed by 
requiring that it vanishes in the vacuum. For the description of hybrid star matter, the results from 
the NIL model have to be merged with a low-density nuclear equation of state. We demand that the 
pressure constant is fixed such that the chiral phase transition coincides with the transition from the 
hadronic model to the NJL model description. Numerically one finds, that the two different pressure 
constants differ only slightly from each other. In the former case one finds a phase transition 
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directly to CFL quark matter, in the latter case two phase transitions appear, first to the 2SC phase 
then to the CFL phase. For a phase transition directly to the CFL phase, the solution is first unstable 
but turns then into a stable one [11]. The new stable solution is another example of the third family 
of compact stars which can appear for a strong first order phase transition [ jlSj [l2| , p^ , [TtI ]. For the 
case of two phase transitions, two kinks appear in the mass-radius curve with stable solutions for 
all configurations. It is interesting that it is possible that there are actually two phase transitions 
present in compact star matter. 

The final state in the evolution of stars with a mass of more than 8 solar masses is a core- 



collapse supernova or a direct collapse to a black hole, see [18] for an overview. The degenerate 
core collapses until normal nuclear matter densities are reached. The repulsion between the nucle- 
ons halts the collapse. A shock front is generated which moves outward but stalls around 100 km. 
The newly bom proto-neutron star emits neutrinos and antineutrinos which deposit energy behind 
the stalled shock. For low-mass progenitors with masses of eight to ten solar masses present su- 
pernova simulations indeed find that this additional energy deposit is enough to lead to a delayed 
explosion mechanism [19]. For more massive progenitor stars, it was recently proposed that the 
standing accretion shock instability (SASI) is crucial for core-collapse supernova explosions after 
about 600ms 

The conditions of core-collapse supernova matter at bounce are energy densities of £ ^ ( 1 — 
1.5)£o> so slightly above normal nuclear matter energy density, temperatures of T ~ 10 — 20 MeV 
and proton fractions of Yp '-^ 0.2 — 0.3 (normal nuclear matter has Yp = 0.5). In the following we 
explore the consequences of a QCD phase transition occurring shortly after the bounce in core- 
collapse supernovae [21]. For simplicity, the MIT bag model is used to describe the high-density 
part of the equation of state, while for the low-density part the supernova equation of state of Shen 
et al. is adopted [22]. The mass-radius relation of neutron stars for an equation of state with a QCD 
phase transition can be drastically different as shown above. The maximum masses are found to 
be 1.56M0 for a bag constant of B^/"^ = 162 MeV and l.SOM© for B^^^ = 165 MeV, so above the 
Hulse-Taylor pulsar mass limit. 

The phase transition line of relevance for astrophysical applications does not coincide with 
the one of ordinary nuclear matter as the conditions are entirely different. The timescales are 
much longer so that strangeness is not a conserved quantum number which makes quark matter 
much more stable. Moreover, supernova matter is neutron-rich so that nuclear matter is unfavoured 
due to the repulsive nature of the asymmetry energy. Both effects lead to a considerably reduced 
critical density for the QCD phase transition, so that the production of quark matter in supernovae 
becomes more likely compared to the situation in heavy-ion physics and can appear already close 
to the conditions at bounce. 

For the situation in heavy-ion collisions, just up- and down-quark matter has to be compared to 
symmetric nuclear matter, the timescales in heavy-ion collisions are just too short to allow for net 
strangeness being produced. Large critical densities are found, several times normal nuclear matter 
densities, so that the production of ud-quark matter is unfavoured for heavy-ion collisions at small 
temperatures and high baryon densities, there is no contradiction with heavy-ion data. To be more 
specific, the freeze-out parameters extracted from statistical models of particle production, /ly „ = 
700 - 800 MeV, Tf,o, = 50-70 MeV for heavy-ion coUisions at SIS energies and /i/.o. ~ 500 MeV, 
Tf o. ~ 120 MeV for AGS energies, are well within the hadronic side of the phase diagram. 
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The equation of state with a QCD phase transition is used as input in a supernova simulation 
which also implements accurate neutrino transport of all neutrino flavours |21]. Shortly after the 
bounce, the QCD phase transition is reached in the core of the proto-neutron star. The formation 
of a core of pure quark matter produces a second shock wave which is moving outwards leading 
to an explosion. Note, that in the standard ID simulation without a phase transition no explosion 
is found, the first shock simply stalls and matter is eventually collapsing to a black hole. The 
supernova simulation runs were performed for lOM© and 15Mq progenitor stars with two different 
bag parameters controlling the onset of the QCD phase transition. We find that the time of the 
appearance of the quark core, the baryonic mass of the compact remnant and the explosion energy 
are significantly sensitive to the location of the QCD phase transition. With the curretn parameter 
sets the quark core appears at tp\, = 200 to 500 ms. If quark matter appears later, the explosion 
energy gets larger, as the density contrast at the accretion shock front increases. Heavy progenitor 
star masses can lead to the formation of a black hole. There are observable implications for the 
QCD phase transition in the neutrino signal of supernovae. When the second shock front produced 
by the phase transition runs over the neutrinosphere a second burst of anti-neutrinos is released. 
The peak location and height is determined by the critical density and strength of the QCD phase 
transition. 

In summary, a first order QCD phase transition at high baryon density can lead to observable 
astrophysical signals involving compact stars and supernovae. Neutron stars with a core of CFL 
quark matter can be stable and can form a third family of compact stars besides white dwarfs 
and ordinary neutron stars. In an exploratory study, we demonstrated that quark matter can be 
formed in supernovae, even shortly after the bounce which produces a second shock with enough 
energy to cause an explosion even in ID simulations. The second shock forms a second peak in 
the (anti-)neutrino signal. Possible implications for the gravitational wave signal and r-process 
nucleosynthesis need to be explored in future work. And much more refined models of QCD 
are needed to explore the consequences of the QCD phase transitions at high baryon densities in 
astrophysics on a more fundamental basis. 

This work is supported by the German Research Foundation (DFG) within the framework 
of the excellence initiative through the Heidelberg Graduate School of Fundamental Physics, the 
Gesellschaft fiir Schwerionenforschung mbH Darmstadt, Germany, the Helmholtz Research School 
for Quark Matter Studies, the Helmholtz Alliance Program of the Helmholtz Association, contract 
HA-216 "Extremes of Density and Temperature: Cosmic Matter in the Laboratory", the Frankfurt 
Institute for Advanced Studies, the Italian National Institute for Nuclear Physics, the Swiss Na- 
tional Science Foundation under the grant numbers PP002- 106627/1 and PP200020-105328/1, and 
the ESF CompStar program. A.M. is supported at the Oak Ridge National Laboratory, managed by 
UT-Battelle, LLC, for the U.S. Department of Energy under contract DE-AC05-00OR22725. 



References 



[1] M. G. Alford, A. Schmitt, K. Rajagopal, and T. Schafer arXiv: 0709 . 4635 [hep-ph 



[2] L. McLerra n and R. D. Pisarski Nucl. Phys. A796 (2007) 83-100, arXiv:0706.2191 
[hep-ph] . 



4 



Astrophysical Implications of the QCD phase transition 



Jiirgen Schaffner-Bielich 



[3] J. M. Weisberg and J. H. Taylor, "The relativistic binary pulsar bl913+16: Thirty years of 

observations and analysis," in Binary Radio Pulsars, F. A. Rasio and I. H. Stairs, eds., vol. 328 of 



Astronomical Society of the Pacific Conference Series, p. 25. 2005. astro-ph/ 04 0714 9 



[4] A. J. Faulkner, M. Kramer, A. G. Lyne, R. N. Manchester, M. A. McLaughlin, I. H. Stairs, G. Hobbs, 
A. Possenti, D. R. Lorimer, N. D'Amico, F. Camilo, and M. Bmgay Astrophys. J. 618 (2005) 

L119-L122, |astro-ph/04117 9 6 . 

[5] D. J. Champion, S. M. Ransom, R Lazarus, F Camilo, C. Bassa, V. M. Kaspi, D. J. Nice, R C. C. 
Freire, I. H. Stairs, J. vanLeeuwen, B. W. Stappers, J. M. Cordes, J. W. T. Hessels, D. R. Lorimer, 
Z. Arzoumanian, D. C. Backer, N. D. R. Bhat, S. Chatterjee, I. Cognard, J. S. Deneva, C.-A. 
Faucher-Giguere, B. M. Gaensler, J. Han, F. A. Jenet, L. Kasian, V. I. Kondratiev, M. Kramer, 
J. Lazio, M. A. McLaughlin, A. Venkataraman, and W. Vlemmings Science 320 (2008) 1309, 



arXiv:0805.2396 [astro-ph; 



[6] talk given at the International Conference on Dense Matter, Hirschegg, Austria, January 19-23, 2009. 
[7] D. J. Nice, I. H. Stairs, and L. E. Kasian A/f Conference Proceedings 983 (2008) 453-458. 



[8] F Ozel Nature 441 (2006) 1 1 15-1 1 17, |astro-ph/0605106 . 

[9] K. J. Pearson, R. Hynes, D. Steeghs, P. Jonker, C. Haswell, A. King, K. O'Brien, G. Nelemans, and 



M. Mendez Astrophys. J. 648 (2006) 1169-1180, 5stro-ph/0605634 



[10] M. Alford, D. Blaschke, A. Drago, T. Kliihn, G. Pagliara, and J. Schaffner-Bielich A^afMre 445 (2006) 



E7-E8, astro-ph/0606524 



[11] G. PagUara an d J. Schaffner-BieHch Phys. Rev. D 77 (2008) 063004, a rXiv :0711.1119 
[astro-ph] . 

[12] B. Kampfer J. Phys. A 14 (1981) L471-L475. 

[13] S. B. Ruster, V. Werth, M. Buballa, 1. A. Shovkovy, and D. H. Rischke Phys. Rev. D 72 (2005) 



034004, hep-ph/050318 4 



[14] D. Blaschke, S. Fredriksson, H. Grigorian, A. M. Oztas, and F. Sandin Phys. Rev. D 72 (2005) 



065020, arXiv:hep-ph/ 05031 94 



[15] U. H. Gerlach Phys. Rev 172 (1968) 1325. 



[16] N. K. Glendenning and C. Kettner Aifran. Astrophys. 353 (2000) L9, astro-ph/9807155 



[17] K. Schertler, C. Greiner, J. Schaffner-Bielich, and M. H. Thoma Nucl. Phys. A677 (2000) 463, 



astro-ph/0001467 



[18] H.-T. Janka, K. Langanke, A. Marek, G. Martmez-Pinedo, and B. Miiller Phys. Rept. 442 (2007) 



38-74, astro-ph/0612072 



[19] F S. Kitaura, H.-T. Janka, and W. Hillebrandt Aifron. Astrophys. 450 (2006) 345-350, 



arXiv: astro-ph/ 05 12 65. 



[20] A. Marek and H. T. Janka arXiv: 0708 . 3372 [astro-ph; 



[21] I. Sagert, M. Hempel, G. Pagliara, J. Schaffner-Bielich, T. Fischer, A. Mezzacappa, F.-K. Thielemann, 



and M. Liehendodei Phys. Rev Lett. 102(2008)081101, arXiv: 0809 . 4225 [astro-ph; 



[22] H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi Nucl. Phys. A637 (1998) 435-450, 



nucl-th/9805035 



5 



